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The plate elements of the NASTRAN code are used to analyze two annular plate 
problems loaded beyond the elastic limit. The first problem is an elastic-plastic 
annular plate loaded externally by two concentrated forces. The second problem is 
stressed radially by uniform internal pressure for which an exact analytical solution 
is available. A comparison of the two approaches together with an assessment of the 
NASTRAN code is given. 


INTRODUCTION 

The piecewise linear analysis option of the NASTRAN code can be used to analyze 
quite complicated elastic-plastic plane-stress problems (refs. 1 and 2). The 
reliability of this code has been well demonstrated in the linear range but not so in 
the nonlinear range of loadings. One major reason is because exact analytical 
solutions for elastic-plastic problems are usually not available for comparison with 
approximate NASTRAN solutions. 

In this paper, the plate elements of the NASTRAN code are ised to solve two 
annular plate problems loaded beyond the elastic limit. The fl st problem considered 
is an elastic-plastic annular plate loaded externally by two concentrated forces. 
There is no analytical solution for this two-dimensional plane-stress problem and the 
NASTRAN code is used to obtain numerical results. The second problem considered is 
an elastic-plastic annular plate radially stressed by uniform internal pressure. 

This problem is chosen because an exact analytical solution is available for 
comparison. For Ideally plastic materials, the stress solution for this statically 
determinate problem was first obtained by Mises (ref. 3) and the corresponding two 
strain solutions were obtained by the present author on the basis of both J 2 
deformation and flow theories (ref. 4). For elastic-plastic strain-hardening 
materials, an exact complete solution was recently reported in reference 5. 

Analytical expressions were derived on the basis of J 2 deformation theory, the Hill's 
yield criterion, and a modified Raraberg-Osgood law. The validity of the above 
solution has been established by satisfying the Budiansky's criterion. 

In the following, the theory of elastic-plastic plate elements as used in 
NASTRAN is briefly reviewed. In its present form, NASTRAN cannot be used for 
problems involving ideally plastic materials. It is shown that this limitation can 
be easily removed by making minor changes. For an elastic-plastic strain-hardening 
material, the NASTRAN solution is reported here and compared with the exact solution. 
The results are presented graphically and an assessment of the NASTRAN code is made. 
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The Cheoretical basis of two-dimensional plastic deformation as used in NASTRAN 
is that developed by Swedlow (ref. 6). In the development, a unique relationship 
between the octahedral stress, t 0 , and the plastic octahedral strain, e Q P, is assumed 
to exist and the use of Ideally plastic materials its excluded. The total strain 
components (e x , e y , e z , and y X y) are composed of the elastic, recoverable 
deformations and the plastic portions (e x p * e y P, c z P » and Y xy P). The rates of 
plastic flow, (e x P, etc.), are Independent of a time scale ana are simply used for 
convenience instead of incremental values. The definitions of the octahedral stress 
and the octahedral plastic strain rate for isotropic materials are: 

t 0 - /(sii z + Zan* + ®22 Z J ' S33 Z )/3 (1) 
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The su is called the deviator of the stress tensor; o x , a y , and T xy are the 
cartesian stresses. The Isotropic material ic assumed to obey the Mises yield 
criterion and the Prandtl fteuss flow rule. The matrix relationship for the plastic 
flow is: 
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where 


2M t (t 0 ) - T 0 /r 0 P 


(5) 


The plastic modulus, Mt(t 0 ), can be related to ti’« slope Ex, of the effective stress- 
strain curve by: 

1 1 1 

3Mf(T 0 ) Ex E 

The total strain increments, obtained by adding the plastic and linear elastic 
parts, are: 


Ue} - ([DPI + [G 1 — 1 ) { Aa } - IGjJ^Uo} 


(7) 
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where [G] is Che normal elastic material matrix and [Gp] Is the equivalent plastic 
material matrix. The matrices [D p ] and [Gp) -1 exist for finite values of or (E?) 
and (Gp) can be obtained numerically. Because this procedure Is chosen In developing 
NASTRAN program, only strain hardening materials can be considered for applications. 
However, It should be noted that even the matrix [Gp] -1 does not exist when Mr or &T 
is equal to zero, the matrix [Gp] may still exist. In fact, the closed form of [Gp] 
has been given In reference 7. We can express this as: 

f s 2 2 2 +2A SYM. 'I 
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If we want to remove the limitation that the use of Ideally plastic materials Is 
excluded, we have to make minor changes in subroutines PSTRM and PKTRM of the NASTRAN 
program . 


TWO-DIMENSIONAL PROGRAM 


Consider a two-dimensional annular plate loaded externally by two concentrated 
forces. Figure 1 shows a finite element representation for one quarter of the 
annular plate. The other part is not needed because of symmetry. There are 198 
grids and 170 quadrilateral elements in this model. The grids (1 through 11) along 
the x-axls are constrained in y-direction while those grids (188 through 198) along 
the y-axis are constrained in the x-dlrectlon. The concentrated force F Is applied 
at the top of the y-axis (grid 198). The thickness of the plate Is 0.1 inch. All 
membrane elements are stress dependent materials. The effective stress-strain curve 
Is defined by: 


where 


o - Ee for c < a 0 
(a/o Q ) n - (E/o 0 ) e for o > a Q 


( 10 ) 


a - (3//2 )t 0 , e - a/E + e p (11) 

n is the strain hardening parameter and the Initial yield surface is defined by the 
ellipse a « a 0 . The Input parameters for the problem are: 

E * 10.5 x 10 6 psl , v ■ 0.3 (Poisson's ratio) , a 0 * 5.5 x 10** psl , n - 9 . 


There Is no analytical solution for the two-dimensional plane-stress problem and the 
NASTRAN code is used to obtain numerical results. First the stress solution in the 
elastic loading range Is obtained and the elastic limit is determined. The corre- 
sponding F* at initial yielding is 753.34 pounds. Then the solution beyond the elas- 
tic limit is obtained in 13 steps with the applied force given by F n m 753,35 (0.95 + 


0.05n) pounds for n - 1,2,... 13. The vertical displacement at the point of applica- 
tion (Vp) and the horizontal displacement (Ujj) at point B (grid 11) are shown in fig- 
ure 2 as functions of the applied force F. The major principal stresses in elements 
near the x and y-axes are shown in figure 3 for F * 600 and 1205.38 pounds, respec- 
tively. The results indicate that the maximum tensile stress occurs at point C. 


ONE-DIMENSIONAL PROBLEM 


Consider a one-dimensional annular plate stressed radially by uniform Internal 
pressure. The plate geometry and material properties are the same as the two- 
dimensional one. Utilizing the condition of axial symmetry, we need only a sector of 
the annular plate for the finite element model. There are only 10 elements with 22 
grids. All grid points are constrained in the tangential direction. The applied 
load is the internal pressure p. If p ■ 1000 psi, the equivalent nodal force at each 
of the two interior grids is Q « 4.62329 pounds in the radial direction. The true 
pressure corresponding to initial yielding for this problem is p Q * 23,571 psi. The 
NAST AN results will depend on the users' choice of element sizes and load Increments. 
In the elastic range, the NASTRAN results based on two finite-element models were 
compared with the exact solution. For a ten— element model, the maximum error is 0.36 
percent in displacements and 0.95 percent in stresses. For a twenty-element model, 
the maximum errors are reduced to 0.09 percent and 0.24 percent, respectively. Since 
we are satisfied with one percent error, the ten-element model is chosen for 
incremental analysis beyond the elastic lir»it. 

In the plastic range, the user has to choose the load increments properly in 
order to obtain good results at reasonable cost. The values of the load factors can 
be normalized so its unit value corresponds to the limit of elastic solution, l.e., 
p 0 » 23,571 psi. The load increments can be uniform or nonuniform. It seems that 
the size of load increments depends on the material properties and sizes of elements. 
In order to determine the influence of load factors on the displacements and stresses 
in the plastic range, four sets of load factors are chosen. The load Increments for 
three of them are uniform with Ap/p 0 ■ 0.20, 0.10, and 0.05, respectively. The 
influence of load factor, p/p Q , on the inside radial displacement, uj, is shown in 
figure 4. The effect of load factors on the major principal stresses in elements is 
shown in figure 5. Vfe also show in these two figures the corresponding analytical 
results. On the basis of these comparisons, we can make the following conclusions. 

In the earlier stages of plastic deformation, larger load increments can be used to 
give very good results. As plastic deformation becomes bigger, smaller load 
increments should be used in order to get a reasonably good answer. However, for 
very large plastic deformation, it seems that we cannot improve the NASTRAN results 
much better by choosing even smaller increments. This is because there are other 
built-in errors in the NASTRAN program, e.g., the linear displacement function is 
assumed. 

If uniform load increments are used, a direct comparison of the analytical and 
NASTRAN results is not available. The solid curves shown in figures 4 and 5 were 
obtained indirectly since the analytical results of the displacements, stresses, and 
pressure were represented as functions of elastic-plastic boundary (ref. 5). Vfe have 
obtained the analytical results when the elastic-plastic boundary is located at a 
radial distance of 0.05, 0.15, 0.25, 0.35, 0.45, 0.55, 0.65, 0.75, 0.85, and 0.95 
inch from the inside surface. The corresponding values of the pressure factor (p/p 0 ) 
are 1.095, 1.266, 1.414, 1.540, 1.646, 1.734, 1 . 805 , 1.860, 1.901, and 1.929, 
respectively. This set of values Is used as load factors in the input deck of the 
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NASTRAN program. Sosa* of the HASTRAN results together tilth the corresponding 
analytical results are shown In figures 6 and 7. A direct comparison of the two 
approaches In the plastic range can be seen. The effect of load factors on the 
distribution of radial displacements Is shown in figure 6. The distributions of 
major principal stresses for three load factors are shown in figure 7. As can be 
seen in figures 6 and 7, a direct comparison of two approaches will support the 
following conclusion. Even if the results in the elastic range are in excellent 
agreement, the differences in the plastic range can be quite big for large values of 
load factors. This suggests sore research efforts should be given to large plastic 
deformation. 


conclusions ORIGINAL PAGE 12 
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Two elastic-plastic annular plate problems have been analyzed by using NASTRAN 
plate elements. One problem is loaded externally by two concentrated forces and the 
other by uniform internal pressure. The NASTRAN results for the second problem have 
been compared with an exact analytical solution. It seems that the NASTRAN code in 
its present form is still a valuable tool because it can be used to solve quite 
complicated plane stress problems provided the plastic deformation Involved is not 
too large. The limitation that the use of ideally plastic materials is excluded can 
be easily removed by making minor changes. 
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Figure 1. Finite Element Model of an Annular Plate. 
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Figure 5. Influence of Load Factor on the Major Principal Stresses 
in Element 1 - Problem 2. 






MAJOR PRIN- STRESS * 10 psi 


ORIGINAL PAGE (9 
OF POOR QUALITY 





Figure 7. Distribution of Major Principal Stresses - Problem 2. 
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